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On Some Expansions for Bessel Integral Functions 

F. Oberhettinger l 

Investigations by van der Pol and Humbert concerning the Bessel integral function of 
order zero are extended to Bessel functions of other kinds and to functions related to Bessel 
functions. 



The Bessel integral function of order zero, 

r* co 
Jto(3)=:-J t- l J {t)dt, 

was introduced by van der Pol [8]. 2 Equivalent [3, 4, 6] with (1) is 

1 f T • 
J%o(x)=- Ci(x sin t)dt. 
ttJo 



(1) 



(2) 



Various properties of this function have been investigated by Humbert [4]. A number of 
definite integrals, especially of the Fourier or Laplace transform type, can be reduced to (1). 
For this reason tables of Ji (x) have been published previously [5]. The numerical computa- 
tion of Jio(x) requires the transformation of the integral in (1) into an expression involving a 
series in either ascending powers (for not too large x) or descending powers (for large x) of the 
variable r. One of these expansions [4, 6] is 



where one can write 



J t0 (tt) =7+ log-+g 2n{n < )2 , 



(3) 



(-D w (*/2) 2 
2n(n\) 2 



= J/~ Vo 



(0-i]*. 



An asymptotic expansion of Ji (x) for large x [5, 7] is listed in (6). The preparation of a 
set of numerical tables by the Computation Laboratory made it necessary to derive expansions 
of types (3) and (6) for functions defined similarly to (1) but with J (t) replaced by another 
kind of Bessel or related function. The established results are listed below. (A list of nota- 
tions and auxiliary formulas is given at the end of this paper.) 



2^(— l) n (x/2) 2n 



2 12 n =i 



2n(n\) 2 



Yio(x)=~j^r l T Q (t)dt=^log |)+T log \ 

[^+l)+^log|] (4) 

*«,)— />«0*-^log |) 2 -7 log H 2 ~£+S SS [^+D+^g 5} (5) 

Both these expressions correspond to (3) and are suitable for not too large x. 
For large x we have the asymptotic expansions of the following functions: 



, Jo(t) 




2 2 2 -4 2 2 -4 2 -6 

•y»2 /yi4 ' /y»6 




1 2 2 2 2 -4 2 2 2 4 2 -6 * , 

X X Z X 5 X 7 



(6) 



1 Guest worker, National Bureau of Standards, from The American University. 

2 Figures in brackets indicate the literature references at the end of this paper. 
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\\- l K,it)dt~-K,{x{ 



j 2 2 2 4,2 2 4 2 -6 
#o(*)| ^+^r+— 3— 



x' a;* 



G + H3 + i 



2 2 2 4 2 2 2 4 2 -6 2 



■]■ 



/» CO 



„tf)-r,(()]«~^S(-i) 



r(2n)!T (2x)- 2 " 
L n! J 2«,+ l' 



£ [ifo(*)-ro(0] <a-?iog (2x)-§y~ -!s(-i) B Si- 



P [I (t)~L (t)] dt—log (2x)~ly 1 2 

Jo 7T 7T T m=1 

Furthermore, we derive the identity 



T—l («,!) 

(2»)!(2n— 1)! 



(«!) 2 



(2x)" 



P r 2 H l (t)dt=(2x)- 1 H 1 (x)+^ r r l H {t)dt- 



(7) 

(8) 

(9) 

(10) 

(ID 



Note that with (6) we know the asymptotic expansion of 



I"'" 1 * 



(0^ for large x by 



virtue of (8). With this and the known asymptotic expansion of H^x) (see (21)), the asymp- 
totic expansion of (11) is likewise known, 
a. Proof of (4) and (5). We write 

/»co n co r* x 

r 1 J v (t)dt= t- l J v {t)dt- t~ l J v {t)dt, 

Jx JO Jo 

and since (using term -by- term integration) 

Jo TO=0 



we obtain by (23) 



m\T(v+m+l)(v+2m) 

(x/2) v ~] ^ (— l) m (x/2)"+ 



X" rw *>*4[ , -ffl-Ss^ 



+ m + l)(*>+2ra) 

This defines an integral function of v. If we now let v tend to zero, using (13) we obtain the 
expansion (3). In order to prove (4), we differentiate the preceeding formula with respect to v. 
If we again let v tend to zero, using (12) and (13) we obtain (4). Finally, the relation 

K,{t)^\iir[J,{it)+iY,{it% 
together with (3) and (4), yields (5). 



b. Proof of (6) and (7). We replace th 
expansion (22), and the results in (6) follow i 


e Lommel functions in (24) by their asymptotic 
immediately. Again, the relation 


K (t)=±i7r[J (it)+iY (it)), 


together with (6), yields (7). (Note that the asymptotic expansion of 




(v/tn 




r» oo 

1 ^ 


Y v {t) 


+ dt 


can be obtained in the same way.) 
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c. Proof of (8). From (21) we have, for v=0 and largo t, 

r 1 [Ho(0-FoW]~-±;(-i)»2- 2 »r^Tr 2 »- 2 . 

7T W =0 L nl J 

Integrating term by term, we obtain (8). 

d. Proof of (9). We use (19) and obtain, interchanging the order of integration, 

P [H (t)-Y (t)]dt=- (~ i-\l+v 2 )- y il~e- xv )dv. 

J0 IT JO 

Because it is necessary to represent the 1. h. s. function of x as an integral of the Laplace 
transform type in order to apply Watson's lemma [1] for the determination of its asymptotic 
expansion, we form, using (14), 

P [H Q (t)-Y Q (t)]dt--log x=- r v- l [(l+v 2 )-X-e- v )dr+- I " e-* v i- l [l-(l + v 2 yy*]dv. 
Jo n njo it J() 

The first, integral on the r. h. s. is a constant equal to 7+log 2. The second integral is of the 
Laplace transform type and can easily be developed asymptotically for large x. We obtain 
first 

v- l [(l+v 2 r H -e- v ]do=\im\ v-\l+tf)-*dv- p- l e-de \ 

JO S-M) LJi J 8 J 

=S ^ j( _ 8) _ Iog[ l±«±^]J 

=7+log 2, by (15). 

The asymptotic expansion of the second integral for large x is as follows. 
We expand 

The application of Watson's lemma [1] gives (9). 

e. Proof of (10). We use (20) to obtain 

P [I (t)L (t)]dt=~ f 1 v-\l-v 2 )- y Kl-e~ vx )dv, 

JO 7T Jo 

and proceed precisely as before, arriving at the intermediate result 

C[Ut)-L Q (t)]dt- 2 log (2x)--y=-f -Ei(-x)+ \\-h-™[\-{l-v 2 yy>]di\- 

JO 7T 7T 7T ^ Jo J 

We expand 



(n\) 2 

The application of Watson's lemma [1] gives (10). (Note that the term —Ei(—x) is 
purely formal and may be omitted, because its asymptotic expansion consists of the product 
of an exponential function with negative exponent and a series in descending powers oix.) 

f. Proof of (11). By means of (18) we obtain 



( a \-*H l (t)dt=-- ( \l-v 2 y> si (px)dv. 

Jx TTJO 
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By partial integration, 

C m l l _ if 1 

t~ 2 H 1 (t)dt = — ~ si (x)+- x l Hi{x)-\— | v * sin * v sin (vx)dv } 

Jx Z Z irj 



but 



f 1 1 1 f °° 

a ! sin x «? sin (vx)dv=^ -k si (x)+- 7r £ l H (t)dt, 



and (11) follows at once. 
Notations: 

«/,(£) =Bessel function of the first kind. 
F„(£)=Bessel function of the second kind. 
I v (t)= modified Bessel function of the first kind. 
K v (t)= modified Bessel function of the second kind. 
Zf„(£)=Struve's function. 
L v (f) = -ie- iv '* /9 H P (it). 
All these functions are Bessel or related functions of the order v. 
T(z) = Gamma function . 

*(s)=r'(z)/r(s). 

7 =-^(l) = 0.577215 (Euler's constant). 

Sn,v(z)=IjommeV8 function. 

P^LH^ [9 - p - 61 ' 62] (12) 

nhr l+yz+ (: 2 -^ ^ (13) 

(Taylor expansion of l/r(l + z) near 2=0). 
Integral representations: 

n oo 

log x= v-\e- v -e~ xv )dv [2, p. 17] (14) 

r* co 

Ei(-<5) = - I v~ l e~ v dv (15) 

/» CO 

Ci(5) = - v~ l cozvdv (16) 

/» CO 

si(5)= — # _1 sin z? dv. (17) 

The last three functions are the exponential integral, the cosine integral, and the sine integral, 
respectively. 

v(v+^H v (t)=2T- l l^hy J\l-^y-^ sm(vt)dv, Tie p>-i t 9 > P- 328 1 (1$ 

H (t)-Y Q (t)=- r e - tv (l+v 2 )- 1!2 dv [9, p. 331] (19) 

I (t)-L (t)=- C e- tv (l-v 2 )- 1/2 dv. [9, p. 338] (20) 

TT JO 
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Asymptotic expansions for large z: 
HM-YJM-*-* cos MQ.y'g /- 1 ^^^^' [9, P- 333] (21) 

S ^,)^-i {i^kiyiL 2 + [(M-l) 2 -^-3)--^] _, , .}. [9> p . 351] (22) 

Integral formulas: 

J"*'- 1 J,(«)<»=2'- 1 r(|+^)/r(l+ir-|^), -Re<Ee M <| [9, p. 391] (23) 



J> 



^=2a:*{ ^_,.,(«)-*^ _ r&-i.o(z). [9, p. 350] (24) 
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